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Analytic expressions for the electromagnetic fields of an ultrashort, tightly 
focused, linearly polarized laser pulse in vacuum are derived from scalar and 
vector potentials, using a small parameter which assumes a small bandwidth 
of the laser pulse. The derived fields are compared with those of the Lax 
series expansion and the complex-source-point approaches and are shown to 
be well-behaved and accurate even in the subcycle pulse regime. We further 
demonstrate that terms stemming from the scalar potential and due to a fast 
varying pulse envelope are non-negligible and may significantly influence laser- 
matter interactions. 
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1. Introduction 

With the rapid development of ultrashort and ultrastrong 
laser pulses^— , laser-matter interactions involving such 
pulses have been the subject of much interest^—. Theo¬ 
retical efforts aimed at describing such pulses have also 
been evolving to meet the need to model the electric and 
magnetic fields of the required pulses^ - — 

In the study of many laser-matter interactions which 
employ a focused laser beam, the laser fields are ade¬ 
quately described by solutions to the wave equation in the 
paraxial approximation^. In the case of a tightly focused 
laser beam, this description is improved upon further us¬ 
ing higher-order corrections to the paraxial solution. For 
example, in the Lax series solution^!! - — , the fields are 
given in ascending powers of the square of the diffraction 
angle e = wo/z r , and provide a good description for a 
focused laser beam of long duration. Here wq and z r are 
the waist radius at focus and the Rayleigh length, respec¬ 
tively. It was further shown that high-order corrections 
to the paraxial solutio n 33 : 34 can have non-negligible con¬ 
tributions for evaluations of direct laser acceleration of 
particles^ - — . 

Another solution of Maxwell’s equations for a focused 
laser beam can be obtained using the complex-source- 
point approach (CSPA)22r—. The fields via CSPA, when 
expanded in terms of the parameter e near the propa¬ 
gation axis of the focused beam, coincide with those via 
the Lax series expansion approximation (LSEA) 41 . For 
a Gaussian beam it was explicitly shown that the sum¬ 
mation of an infinite Lax series reduces exactly to the 
solution of CSPA—. The CSPA solution for a propagat¬ 
ing spherical wave is singular at the point r = ±i(z+iz r ), 
with the coordinate in the propagation direction z and 
the radius r, which limits its utility far from the propaga¬ 
tion axis and causes problems in numerical simulations. 
However, a method for the regularization of CSPA solu¬ 
tions was foun d 42 : 43 . Rather than a propagating spheri¬ 


cal wave, one may use the field of a simple point source 
and a point sink placed at the origin, which constitutes 
a standing spherical wave. Then the focused beam, ac¬ 
cording to CSPA, is free from singularitie s 24 : 42- — . The 
latter field was shown to form when using a laser beam 
focusing with a parabolic mirror. 

However, the use of LSEA, or analogously CSPA, is 
limited in that the pulse temporal shape is factorized in 
the field expressions and, consequently, it does not take 
into account distortions of the pulse during propagation. 
While this may work well for rather long multi-cycle laser 
pulses, in which the change in pulse-shape is not appre¬ 
ciable, yet for many important current and future devel¬ 
opments involving ultrashort and tightly-focused pulses 
of laser light, the space-time evolution and distortion 
of the pulse-shape ought to be taken into account^ - —. 
An exact analytical solution of Maxwell’s equations is 
known which describes a focused laser beam of a spe¬ 
cial form with a cosine-squared pulse shape of arbitrary 
length— However, this solution in the form of a double 
sum is inconvenient for analytical treatments. The 3D 
vectorial Maxwell’s equations are solved in22 for a laser 
pulse, yet, only in the paraxial approximation. The ex¬ 
tension of CSPA to describe pulsed laser light has also 
been implemente d 15 ' 24 ’ 28 ' 47 . Nonetheless, analytical re¬ 
sults have been obtained mostly for isodiffracting laser 
pulse s 15 ' 27 . 

The LSEA approach was generalized to describe the 
spatio-temporal evolution of an ultrashort and tightly- 
focused laser beam by solving the wave equation with 
the use of two small parameters, namely, the spatial 
diffraction angle 1/koWo and the temporal spreading pa¬ 
rameter I/co’oTq, where tq is the laser pulse duration, to o 
the frequency, and ko the wave number (coo = cko in 
vacuum}^ - —. The corrections to the paraxial and slowly 
varying envelope approximations for the fields of the fo¬ 
cused radially polarized laser pulse of the TMoi Gaussian 
mode were derived in? 2 with a purely time-domain ap- 
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proach. It was shown that due to those corrections the 
longitudinal field decreases quickly along the propagation 
direction. 

About two decades ago, Esarey et. proposed a way 
to describe the spatio-temporal evolution of an ultrashort 
and tightly-focused laser beam by solving the wave equa¬ 
tion with the use of a single small parameter, namely, the 
relative bandwidth of the laser pulse vj = Aw/wo -C 1. 
This solution describes the field of an ultrashort laser 
pulse when the duration of the pulse fulfils the condition 
r 0 < z r j c. The advantage here is that the ultrashort 
and tightly-focused laser beam is described by a com¬ 
pact expression in terms of an expansion over a single 
small parameter, in contrast to the approach 0 fi£r—i£S. 
which uses expansions over two small parameters. In^, 
the laser pulse propagating in a plasma was considered 
and the transverse fields (with respect to the propaga¬ 
tion direction) were calculated from a transverse vector 
potential alone. However, the description of the longitu¬ 
dinal components of the laser field, which become very 
significant for tightly focused laser beam s 22 ! 30 , requires 
the introduction of a scalar potential. Moreover, when 
the laser field is described by scalar and vector poten¬ 
tials, the scalar potential will alter the transverse fields, 
as well. 

In this paper we derive analytic expressions for the 
electromagnetic fields in vacuum of an ultrashort and 
tightly focused propagating laser pulse of linear polar¬ 
ization and of a Poisson-like spectrum. The fields are de¬ 
rived using the approach o^, but combining both contri¬ 
butions of vector and scalar potentials to describe prop¬ 
erly the longitudinal field of the ultrashort and tightly fo¬ 
cused laser beam. The field expressions include the first 
order corrections with respect to the small bandwidth 
parameter vo. 

The structure of the paper is as follows. In Sec. [2] 
the calculation of the fields is outlined. Comparison of 
the obtained field structure with the CSPA and LSEA 
solutions for the pulses are discussed in Sec. [3] And, our 
conclusions will be given in Sec. [Tj 

2. Calculation of the fields 

Expressions for the electromagnetic fields propagating in 
vacuum may be found from vector and scalar potentials 
A and $, respectively, which satisfy the equations 


and currents^. Note that, assuming 

$ = (j) 0 (t>{x , y , 0 , t) exp [i(k Q z - w 0 t)], (4) 


where cf> 0 is a constant, and ko = 2 tt/X 0 a central wave 
number, corresponding to the wavelength Ao, the Lorentz 
condition yields 


c 2 V ■ A 

iojQ — (d(j>/dt)/(j>' 


( 5 ) 


In the special case of a vector potential polarized linearly, 
say along the cc-direction, the vector and scalar potentials 
satisfy identical (scalar) wave equations. As such, the 
potentials can only differ by multiplicative constants. We 
may, therefore, write 


A = xaoa(x, y,z,t) exp[z(/coz — wot)], (6) 


where x is a unit vector in the polarization direction and 
ao a constant. Now, setting <f> = a the scalar potential is 
turned into 


c 2 V ■ A 

iu>o — (da/dt)/a 


( 7 ) 


Thus, the problem of finding the scalar and vector poten¬ 
tials reduces to finding a solution for a(x,y,z,t)- Then, 
expressions for the fields E and B may, respectively, be 
found from 

BA 

E = —— - V<f>, and B = V x A. (8) 
dt w 

In the absence of a scalar potential function $, the elec¬ 
tric field E will have the same polarization as the vector 
potential A , whereas the magnetic field B will not be 
altered by introduction of the scalar potential. 

Following Esarey et. a&, a change of variables is first 
applied to £ = z — ct, y = (z + ct)/ 2. Consequently, the 
wave equation, Eq. (1), becomes 



A = 0. 


Inserting Eq. JS]) into Eq. m , it reads 


( 9 ) 



2 ikr 


drj 


d 2 

+ BC,dy 


a (r, £, V) = 0 , 


( 10 ) 


1 d 2 A 



V 2 A 

c 2 dt 2 

= 0 , 

(1) 

.. 1 9 2 <h 



V 2 $ 

c 2 dt 2 

= 0 , 

(2) 


provided the Lorentz condition (SI units are used), 

1 

v - 4 + ?ar = °- < 3 > 

is satisfied, simultaneously. They are equivalent to the 
full set of Maxwell’s equations in the absence of charges 


where r = ^x 2 + y 2 . This equation is exact, i.e., no ap¬ 
proximations have been assumed yet. When the last term 
is neglected in Eq. m, it becomes similar to the parax¬ 
ial equation^, which has an analytical solution. Using 
the scaled variables r = r/wo, fj = T]/z r , and £ = C/ T 0 i 
Eq. ifTOl) reads 

(vi+«A +te A_) o ( f ,<,, ; ) =0 . (ii) 

The last term is proportional to a small parameter vj = 
(wqTo)^ 1 . Taking into account this term perturbatively, 
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the final solution can be derived as a series with powers 
of vo. Most easily the latter can be carried out by Fourier 
transforming Eq. (1101) with regard to fy which yields 


+ 2 i (fc 0 + k) — 


ak (r, k, rj) = 0 , ( 12 ) 


where 


1 

a k {r,k,ri) = —= a (r, C, rj) exp(-ik()d( (13) 

v j — oo 


is the Fourier transformation of a (r, rj). Equation (1121) 
has the same mathematical structure as that of the 
paraxial wave equation, which has explicit solutions for 
different modes. For instance, taking the lowest-order 
Gaussian mode, the Fourier transform of the laser enve¬ 
lope reads' 


a fc (r, k, rj) = fk exp ip k (r, k 0 + k, rj ), (14) 


with 

1 P 2 

ip k (r,k 0 + k,r]) = --ln(l + a 2 k ) - ——— 

2 1 + ia k 

—itan _1 afc, (15) 

where a k = v/ z rk, z rk = (fco + k)w 2 /2 , p = r/w 0 , w 0 
is the laser focal radius, and fk the initial axial envelope 
profile with the focus taken at 77 = 0. 

The laser pulse envelope a (r, rj) is calculated as the 
inverse transform of Eq. 


The inverse Fourier transform can be carried out an¬ 
alytically if , 0 fc(r, fco + fc,p), viewed as a function of 
k' = k + ko, is Taylor-series expanded around the central 
wavenumber kg, as in 


fk exp(^) exp(ikj)dk. (16) 


j>k(r, k 0 + k, rj) 

V4 n) (r, ko,rj) 


00 , LM 

" (a k 0 ,T])—, 

n—0 n - 

d n ipo(r, k 0 ,T]) 

dkjf 


(17) 

(18) 


where rpo(r,ko,rf) = tpk{r,k 0 + k,rj) |fc=o- For instance, 
the first two terms in the series expansion above are 


( 0 ) 

1 

ln(l + a 2 ) 

P 2 

— i tan 1 a, 

(19) 

0 

2 

1 + ia 

( 1 ) 

ia 

1 

P 2 


( 20 ) 

0 — 

k 0 

1 + ia 

(1 + ia) 2 

1 


where a = t]/z r , and z r is the Rayleigh length of the laser 
beam. 

The expansion of Eq. da is valid as long as k/ko -C 
1. From the uncertainty relation (Ak) (At)) ~ 1, with 
Afc ~ k and A£ ~ cto, one obtains k ~ 1/ctq. There¬ 
fore, the fields derived from this approach are valid when 


To l/u>o, which is fulfilled for a pulse length of at least 
half a cycle. Thus, Eq. ED represents the sought series 
expansion over the small parameter vo. 

The expansion Eq. ED is applicable in the space-time 
region where a < 1. Taking into account that a = 77 / z r = 
(z + ct )/2 = z/z r — £/( 2 z r ), the range of applicability of 
the obtained solution may be estimated by z < z r and 
cto < z r , while the latter reads, equivalently, 

w 0 t 0 < (k 0 w 0 ) 2 . (21) 

The laser pulse envelope a(r, £, rj) can, in principle, be 
obtained by an inverse Fourier transform, i.e., by insert¬ 
ing Eq. ED into Eq. ED- Furthermore, the fields E and 
B may then be calculated via Eqs. © and ©. 

A. Laser pulse with Gaussian spectrum 


For calculating the fields of an ultrashort tightly focused 
laser pulse one needs to specify the spectral distribution 
of the field fk in Eq. ED- Simple analytic expressions 
will be obtained when the spectrum of the laser pulse 
is a Gaussian. Although, in this case it is not possible 
to avoid unphysical negative and zero frequency compo¬ 
nents, it still can be applied when the pulse’s bandwidth 
is narrow enough, i.e., Aw -C uiq. More complex solutions 
based on Poissonian spectra are free from the problem 
above and discussed in the next section. 

For the lowest order solution, 

1 r°° 

a(0) (r, C, rj) = exp(v4 0) )^= / f k exp(ik()dk 

V J —00 

= exp(4 0) )/(O, (22) 

since is independent of k. The laser envelope /(£) 
is the inverse Fourier transform of f k ■ The solution of 
Eq. ED i s singular at k = — ko- To avoid the physically 
uninteresting situation arising from k = —ko, when the 
laser field is axially uniform, see Eq. ED, the following 
Gaussian solution of fj is considered^ 



(23) 


where, a = 2 ^hi 2 ’ an< ^ ^ = cr ° e< f ua ^ s full-width-at- 
half-maximum (FWHM) of the laser pulse. The inverse 
Fourier transform of fk gives 



(24) 


With Eqs. ©, (1221) and (1241) . the vector potential of 
the zeroth-order Gaussian solution reads 


A = xao/3 I 1 + 




exp (ikoC — 


k 2 r 2 


,(25) 


where we have used the ansatz 


exp(^ 0) ) 

P 


P exp (- Pp 2 ) , 
1 

1 + ia 


(26) 

( 27 ) 
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Following are the electric and magnetic field components, 
calculated via Eq. ([S]) from Eq. (B5l) 


-Err = 


Ey - 


E z = 


Ei 
Ci 
Eifxy 


Cizl 
Ei fix 

~C^ r 

+ 


-C? + 
p + 


/3 2 x 2 


ipx 2 


Z 2 k 0 Zr CiZr (C 2 - C/2) 
iz r 


Zf ( 2 

IZ. 


Ci{z+iZ r -f/2) 2 _ 
ip 2 r 2 ( k 0 C\ if 

k 0 


4 z 2 a 2 

C 2 z r 


(a 2 - ik 0 f) Ci 


B t = 0, 


B v — Ei ( —i + 


if iP 2 r 2 k 0 C 


1 


2kozr 4 z 2 a 2 ia 2 + kof 


B z = 


Eipy 


(28) 


where 


Ei 


Ci 


C 2 

C 3 


EoP f, , ik 0 f\ ( i a i • 

+ ~^r) ex P ^oC + *v>o - 

. ikjr 2 - C 3 _ 1 fcoC 

C 2 kof + ia 2 a 2 ’ 

C 3 - 2ik 2 r 2 1 1 

Ci + (kof + ia 2 ) 2 + a 2 
2k 0 (z + iz r — C/2), 


Me! 

2a 2 



(29) 


y>o is the carrier envelope phase and Eq = koao the 
field amplitude. The various components are scaled to 
make E/Eq = 1, at x = y = z = t = 0, employ¬ 
ing a normalization factor Fn = i/ ( koz r Co ) — ico, with 
co = 1 + 1 / ( 2 k 0 z r ) + l/a 2 . 

The fields calculated from a vector potential alone^ are 
clearly transverse, with 


E x = -EiCi, 


(30) 


and E y = E z = 0. The magnetic field components, 
however, would be the same as in Eqs. (12811 . Thus we 
see that the scalar potential, which is necessary for de¬ 
scribing the longitudinal fields, also modifies the trans¬ 
verse fields and, in this case, the normalization factor is 
F n = i[l + 1/(2 kzr) + 1/ a 2 ]. 

B. Zeroth-order Poisson solution 


For ultrashort laser pulses the Gaussian spectrum fk con¬ 
tains unphysical negative frequencies. A spectrum which 
does not contain negative frequencies has to be chosen to 
adequately describe ultrashort pulses. Here, we choose 
the Poissonia n 49 ! 50 


fk 



g—sk/ko 
r 0 + 1 ) 


H(k), 


(31) 


where s is a real positive parameter, T(x) the gamma 
function, and H(k ) the unit step function which ensures 


that the pulse does not exhibit negative frequencies. Note 
that the propagation term e 1 ^ is included therein. The 
parameter s determines the laser pulse duration tq (the 
root mean square of the temporal distribution |/(t)| 2 ) 


to 


2s 

ojqi/ 2 s — 1 


(32) 


Taking the inverse Fourier transform of fk given by 
Eq. (1311) . the laser envelope is obtained as 


/ (C)=( 1 -!M) S \ ( 33 ) 


With this, the zeroth-order vector potential takes the 
form 


A (°) 


xao/3 




— s — 1 

exp(-/V). 


(34) 


Finally, the following field components may be obtained 
„ _ \ (2fx 2 i p 2 x 2 \ 

E x Ei - 7-v 9 ( 2 ^ n 2 ) 1 D 1 

iL> iWq \ Wq DiWqZ t J 


Ey 


EiP 2 xy / _ 2/3 \ 


iDiu>Q 

Ei2/3x 

Diwl 


E z = 

By = Ei ( iD 
B z = 
where 


-Di - 


DiZ r J 
2k 0 s' 


D, 


s — ikof Di 
2ikos' 


s — ikof 


2Eif3y 


(35) 


Ei = 


D 2 = - 


EoP I 

\ s 1 

F n ' 

V s J 

P 

P 2 P 2 k 0 s 1 

2 z r 

2 z r s — iko 1 


exp {~Pp 2 + iipo), 


P 2 , P 2 P 


, 2«3 


k 2 s’ 


4 zf 2zf (g — ikof) 


2 ’ 


(36) 

(37) 


and, s' = 1 + s. Also, the normalization factor Fn = 
—do — 2 /dou>o, and do = i/(2z r ) — iko (1 + s) /s. 

C. First-order Poisson solution 

The first-order solution of a(r, f, 77 ), using Eqs. (TTBl) and 
(ED. can be written as 


a ( r , C, v) = 


1 


y/2n J - c 
1 r 


fk exp(^Q 0) + kifl) exp(ik()dk 
fk exp(t/>Q 0) ) exp[*fc(£ - iif { o ] )\dk 


y/2n J- 

= /(C)(C - ii>o ) ) ex p(4 0) )- 


( 38 ) 
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Accordingly, the first-order Poisson-based vector poten¬ 
tial A reads 


A 1 - 1 ) = xao/3 


ik 0 (c~ii/j Q 1] ) 


exp(-/3 P 2 ), 


(39) 

with -0Q 1 ' 1 given by Eq. (15U1) . Using Eq. (1551) . the electro¬ 
magnetic fields are 


E x = 


iE\ 

2 Zf 

+H 4 X 


, , His'\ 2 Et/3 


1 + 


i/3s'a 

H 2 s 


Ey = H A y , 

Eix/3 


E, = 


o-a • 2 Q 2 /3s (z + 3/3a — p 2 P 2 a) 

- ip p - - - 

H 2 s 


H 3 wlz r 

is'(4:kz r — Hi) | /3(2 + s)s'(4fcz r — Hi)a 

Ihs + 

E lX /3 3 * 


2 H 2 z?w 2 


1 + 


i/3s'a 

H 2 s 


H 2 s 2 

s' (2/3 — H{)(Akz r — Hi) 
Fffs 2 /3 2 


+1 - V/3 ~ 2s '( 1 ~ 2 P 2 P ~ + 3fp 2 /3 2 a) 


ff 2 s 


By — 


iEip 

2z r 


2o 1 , s'(4fcz r - Zfi) 

P P- H- 77-5 - 

H 2 ps 


ZEiPy L + */3sAA 
w o V # 2 s y ’ 

'"l - ik 0 (( + /3 


B 
where 
Ei = 


(40) 


n — s —1 


F n l 
exp(-/3p 2 + iipo), 

Hi = 2 kz r — /3 + p 2 /3 2 — 2ip 2 /3 3 a + */3 2 a, 

H 2 = 1 — (c + # 0 ^) / s > 


H 3 = 


da/dct 
a 

W f 2n_ U _ iH l S ' 

2z r {P P J 2H 2 z r s' 


(41) 

(42) 

(43) 


(44) 


Hi = 


4:Ei/3 2 x 

H 3 w$ 


1 | 2i(3s'a 2 f3 2 s'a(2 + s) 

1 T 77 77o t 


ff 2 a 


Hn S 2 


/3(/Uf s 2 — His'a — iH 2 s's — 2H 2 s's/3a) 

2H%H 3 z r s 3 X 

(i? 2 s + is'/3a)], (45) 

s' = s + 1, and the normalization constant is 


F n = 


2 Zr 


^ s'(2kz r — 1) 


iiZrS _. (46) 

Wq(1 — 2 kz r s’) 


3. Discussion 

Typical distributions of the electric field components 

are shown in Fig. [T] Figs. Q] (a)-(c) are based on 


Eqs. (1281) . derived from combined vector and scalar po¬ 
tentials (VSP) while Fig. []Jd) is based on Eq. (15D1) or 
fields derived from a vector potential (VP) alone. Note 
that Figs. [T] (a) and (d) are hardly distinguishable for the 
parameter set used, on account of the fact that the ad¬ 
ditional terms, except the first one in the E x component 
of Eq. (051) . constitute a very small correction to the first 
term, i.e., the right side of Eq. I®. Easily visible devia¬ 
tions may be reached when subcyle pulses are employed 
instead. On the other hand, E z and E y , illustrated by 
Figs. [T] (b) and (c) are solely due to the VSP approach 
and are absent from the VP analysis. Those components, 
although more than two orders of magnitude weaker than 
their associated E x component, yet play a significant role 
in the treatment of interaction of matter with such a laser 
pulse in vacuum. 

Figure [5] compares and contrasts our Gaussian solu¬ 
tions for ultrashort focused laser pulses, see Eq. (l28l) . 
with those via CSPA3 and LSEA2&. In Fig. |U (a) 
the pulse duration is relatively long L = 5Ao, while in 
Fig. [5] (b) the pulse is of subcyle duration. In both cases 
the laser beam is tightly focused with a waist size of 
u>o = Ao- It is hard to distinguish the fields stemming 
from all three approaches for the case of a long laser 
pulse {L = 5Ao). The slight difference between the CSPA 
and LSEA solutions, on the one hand, and ours, on the 
other, in this case is due to the fact that the applica¬ 
bility condition Eq. m of our solution is at its limit, 
namely, uiT 0 /(kw o) 2 ~ 0.8. In the case of a subcycle 
pulse ( L = Ao/2), the variations of E x at focus with 
time based on the CSPA and our Gaussian solution are 
almost identical, whereas those based on the LSEA are 
clearly different due to temporal diffraction effects (the 
LSEA solution is valid only for pulses of infinitely long 



x/\q zj Ao 


Fig. 1. (Color online) Contour plots of electric held com¬ 
ponents normalized by Eq at t = 0. (a)-(c) are based 
on the Gaussian held of Eqs. (l28l) . (a) and (b) show 

E x and E z in the xz —plane (y = 0) and (c) shows E y 
in the xy —plane {z = 0). (d) Same as (a) but based 
on the Gaussian held, which is derived solely from the 
vector potential, see Eq. ®. The parameters used are 
wq = L = Ao = 1 pm. 
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ci/cr 


Fig. 2. (Color online) The time variation of the normal¬ 
ized transverse laser field E x at focus. The pulse lengths 
are (a) L = 5Ao and (b) L = Ao/2. The other parame¬ 
ters are the same as in Fig. [I] The blue dashed curves 
show our solution of Eq. (|28|) . while the solid red and 
dot-dashed black curves are based on fields of CSPA and 
LSEA, respectively. Please note that: (a) all three curves 
are close to each other with merely slight differences; (b) 
the blue dashed and solid red curves are overlapping. 


duration). 



Fig. 3. (Color online) Time variation of the normalized 
E x at focus. The pulse lengths are (a) L « 0.6Ao (s = 2) 
and (b) L ss 2.7Ao (s = 50) with wq = 5Ao, while the 
remaining parameters are the same as in Fig. [T| The 
red solid and blue dashed curves present the fields of 
zeroth-order Gaussian and zeroth-order Poisson solution, 
respectively. 

In Fig. [3l we compare the fields of pulses with differ¬ 
ent spectral distributions, namely, Gaussian and Poisso- 
nian. Note that the red solid and blue dashed curves 
present the fields of zeroth-order Gaussian and zeroth- 
order Poisson solution, respectively. For a subcycle laser 
pulse of L « 0.6Ao (s = 2), see Fig. EJa), the field 
with Gaussian spectrum is significantly different from the 
Poisson solution. Therefore, the Gaussian spectral pro¬ 
file for the pulse is not suitable anymore, and the Poisson 
field should be employed to simulate ultrashort pulses. In 
Fig. [3](b), for a several-cycle laser pulse with L ss 2.7Ao 
(s = 50), the Gaussian’s field overlaps with the Pois¬ 
son solution near the pulse center, and they very slightly 
differ from each other around the pulse’s tail ct/a > 8. 


Therefore, the Gaussian spectrum is assumed to be an 
acceptable solution for this case. 



Fig. 4. (Color online) Time variation of the normalized 
E x at focus, (a) Wq = Ao, and (b) Wq = 0.5Ao with 
L « 0.6Ao (s = 2), while the remaining parameters are 
the same as in Fig. |1] The blue solid and red dotted 
curves present the fields of zeroth- and first-order Poisson 
solutions, respectively. 


We compare the zeroth- and first-order Poisson solu¬ 
tions for a subcycle laser pulse with different laser focus¬ 
ing sizes in Fig. [I] Note that the blue solid and red dotted 
curves present the fields of zeroth- and first-order Poisson 
solutions, respectively. For a tightly focused laser pulse 
of wo = Ao, see Fig. Oja), the first-order solution of the 
Poisson field overlaps with the zeroth-order solution very 
well. Therefore, the zeroth-order Poisson field suffices to 
simulate such pulses. In Fig. [3Kb), for an ultra-tightly fo¬ 
cused laser pulse with wq = 0.5Ao, the first-order solution 
of the Poisson field overlaps with the zeroth-order solu¬ 
tion near the pulse center as well, however, they are dis- 
tinguishably different from each other around the pulse’s 
tail ct/a > 1. Consequently, in this case the zeroth-order 
solution is not sufficient anymore, and the first-order cor¬ 
rection should be taken into account. 

4. Conclusion 

Following the approach of2, we have derived the field ex¬ 
pressions of an ultrashort tightly focused laser pulse prop¬ 
agating in vacuum, using an expansion over a single small 
parameter which assumes that the relative bandwidth of 
the laser pulse is not large. The basic solution is derived 
using a Poissonian initial spectral distribution for the 
pulse, and it is applicable for pulse lengths not larger than 
the Rayleigh length. The accuracy of the solution is ver¬ 
ified in comparison with that of CSPA of2£. It is demon¬ 
strated, with the introduction of a scalar potential, that 
the fields do not only develop new axial components, but 
also modify the transverse ones quite significantly. More¬ 
over, deviations of the derived fields from those based 
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on the LSEA are shown to be non-negligible when the 
pulse duration is less than a few cycles, and substantial 
for subcycle pulses. It has been shown that the Gaussian 
spectral distribution cannot be applied in the case of sub¬ 
cycle laser pulses, because they yield pulses differing sig¬ 
nificantly from those via the Poissonion-spectrum-based 
ones, with the same bandwidth. Compared with solu¬ 
tions of the wave equation expanded in power-series of 
two parameters, the diffraction parameter ( 1/kwo ) and 
the time dilation parameter (1 /wtq), respectively, our so¬ 
lution is more compact, cf., e.g.ji^, and may be of utility 
for current and future applications in laser acceleration 
and other ultrafast physics. 
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